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Pheochromocytoma (PCC) and abdominal paraganglioma (PGL) are neuroendocrine tumors that present with clinical 
symptoms related to increased catecholamine levels. About a third of the cases are associated with constitutional muta- 
tions in pre-disposing genes, of which some may also be somatically mutated in sporadic cases. However, little is known 
about inactivating epigenetic events through promoter methylation in these very genes. Using bisulphite pyrosequenc- 
ing we assessed the methylation density of 11 PCC/PGL disease genes in 96 tumors (83 PCCs and 13 PGLs) and 34 normal 
adrenal references. Gene expression levels were determined by quantitative RT-PCR. Both tumors and normal adrenal 
samples exhibited low methylation index (IVletl) in the ECLNl (PDH2j, MAX, MENl, NFl, SDHB, SDHC, SDHD, SDHAF2 (SDH5), 
and 77MEM727 promoters, not exceeding 10% in any of the samples investigated. Aberrant REfpromoter methylation was 
observed in two cases only. For the VHL gene we found increased IVletl in tumors as compared with normal adrenals (57% 
vs. 27%; P < 0.001), in malignant vs. benign tumors (63% vs. 55%; P < 0.05), and in PGL vs. PCC (66% vs. 55%; P < 0.0005). 
Decreased expression of the VHL gene was observed in all tumors compared with normal adrenals (P < 0.001). VHL Metl 
and gene expressions were inversely correlated (R = -0.359, P < 0.0001). Our results show that the VHL gene promoter 
has increased methylation compared with normal adrenals (IVletl > 50%) in approximately 75% of PCCs and PGLs investi- 
gated, highlighting the role of VHL in the development of these tumors. 



Introduction 

Pheochromocytoma (PCC) is a neuroendocrine tumor 
located in the adrenal medulla. Abdominal paraganglioma 
(PGL) is a highly related tumor occurring outside the adrenal 
medulla, most commonly within the sympathetic paravertebral 
ganglia. 

The clinical profile of PCC and PGL is often characterized 
by increased catecholamine secretion and associated symptoms 
— most pronounced in tumors with epinephrine overexpression.' 



The majority of patients present with common endocrine-related 
manifestations such as headache, palpitations, tachycardia and 
excessive perspiration,^ but asymptomatic tumors are also com- 
mon and many are found only at autopsy.^ 

The last few decades have witnessed major discoveries relat- 
ing to the genetic background of these tumors. Approximately 
one third of PCC and PGL patients carry a constitutional muta- 
tion in a predisposing gene including EGLNl (PDH2), EPASl 
(HIF2A), KIFIB^, MAX, NFl, RET, SDHA, SDHAF2 (SDH5), 
SDHB, SDHC, SDHD, TMEM127, and Vm."' Somatic 
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Table 1. Results of methylation analyses in the tumors and normal adrenals 



Gene 
symbol 


Gene name 


Tumors 


Adrenals 


Tumors vs. 
Adrenals 
P value 


n 


Metl 
mean % 


Metl 
range % 


n 


Metl 
mean % 


Metl 
range % 


EGLN1 


Egl nine homolog 1 


95 


1 .6 


0-5.6 


1 2 


1 .2 


0-2.4 


n.s. 


MAX 


Myc associated factor X 


96 


1 .0 


0-4.7 


6 


0.9 


0-3.2 


n.s. 


MEN1 


iVlultiple endocrine neoplasia type 1 


95 


1 .6 


0-5.7 


1 1 


1 .6 


0.9-4.6 


n.s. 


NFl 


Neurofibromin 1 


96 


1 .0 


0-3.8 


4 


0.6 


0-1 .5 


n.s. 


RET 


Rearranged during transfection proto-oncogene 


96 


2 8 


0-28 8 


1 2 


2 9 


1 Q_3 g 


n.s. 


SDHAF2 


Succinate dehydrogenase complex assembly factor 2 


96 


0.8 


0-4.0 


12 


0.8 


0-2.5 


n.s. 


SDHB 


Succinate dehydrogenase complex, subunit B 


93 


1.1 


0-1.7 


15 


1.5 


0.6-2.0 


n.s. 


SDHC 


Succinate dehydrogenase complex, subunit C 


95 


0.3 


0-1.5 


18 


0.6 


0-2.3 


n.s. 


SDHD 


Succinate dehydrogenase complex, subunit D 


96 


1.5 


0-3.5 


15 


0.7 


0-1.5 


n.s. 


TMEM127 


Transmembrane protein 127 


96 


1.5 


0-9.0 


11 


2.3 


0.9-5.0 


n.s. 


VHL 


Von Hippel-Lindau tumor suppressor 


96 


56.5 


19.1-88.6 


34 


27.0 


14.4-43.8 


P< 0.0001 



n.s., not statistically significant; IVIetl, methylation index 



mutations in some of these genes also occur, most commonly 
involving 7VF7.'°'" 

Studies of DNA methylation have revealed CpG hypermethyl- 
ator phenotypes in PGLs carrying a SDHB,^^'^^ or SDHx mutation 
as was also demonstrated in a mouse model.''' However, possible 
somatic inactivation of the PCC/PGL predisposing genes by pro- 
moter methylation has so far not been fully determined. Gascon 
et al. reported unmethylated promoters for SDHB, SDHD, and 
VHL in 35 PGGs," and Dammann et al. described unmethyl- 
ated VHL promoter in 25 PGGs.'* Later, Welander et al. reported 
unmethylated promoters for KIFIB^, MAX, NFl, RET, SDHA, 
SDHAF2, SDHB, SDHC, SDHD, TMEM127, and VHL in 42 
sporadic PCGs." These three studies have employed methyla- 
tion-specific PGR (MSP) for qualitative identification of hyper- 
methylation. To unambiguously evaluate the gene promoter 
methylation status in these tumors, we employed quantitative 
Pyrosequencing in a large panel of PGGs and PGLs. Findings 
were subsequently compared with gene expression levels and eval- 
uated against tumor and patient phenotypes. 

Results 

Increased promoter methylation density of VHL in PCC 
and PGL 

CpG methylation density was quantified using Pyrosequencing 
for the promoter regions of 11 genes implicated in PGC/PGL 
{EGLNl [PDH2], MAX, MENl, NFl, RET, SDHAF2 [SDH5], 
SDHB, SDHC, SDHD, TMEM127, and VHL) in 96 tumors 
(83 PGGs, 13 PGLs) and 34 normal adrenal samples. For each 
sample and gene a methylation index (Metl) was defined as the 
mean methylation of all GpG dinucleotides covered. Overall, the 
methylation assessment for all genes except VHL showed no or 
very little methylation. With regards to EGLNl, MAX, MENl, 
NFl, SDHAF2, SDHB, SDHC, SDHD, and TMEM127, Metl 
values below 10% were obtained for all tumors and normal 



adrenals (Table 1). One single GpG site (case PGG-56) within 
the TMEM127 gene showed an increased methylation level of 
15%. For the RET gene increased methylation was recorded in 
two PGLs (PGL-102 and PGL-113) with Metl values of 17% and 
29%, respectively, while all other tumors and normal adrenals 
showed Metl below 10%. 

By contrast, increased methylation was frequently observed 
for the VHL promoter (Fig. SI; Table 2). Gomparing VHL meth- 
ylation between sample groups (Fig. 1) revealed increased Metl 
in tumors as compared with normal adrenals (mean Metl 57% 
vs. 27%; P < 0.0001). DNA from 4 normal adrenal medullary 
samples (NAM-1— 4) also showed low methylation levels (mean 
Metl = 31%) comparable to the normal adrenals. Altogether 
72/96 (75%) of the tumors had a Met! > 50% while the Met! 
of the normal adrenals did not exceed 44% in any single case 
(Table 1). We further observed a difference in methylation level 
between specific GpG sites within the analyzed VHL promoter 
sequence (Table 2; Fig. 2; Fig. SI). This was particularly pro- 
nounced at GpG 1 and 2 (Ghr3: 10 183 196 and Ghr3: 10 183 198 
respectively). Specifically, we observed a statistically significant 
(T < 0.0001, Table 2) higher methylation level in tumors as com- 
pared with normal adrenals at GpG 1 (mean methylation 80% 
vs. 37%), GpG 2 (72% vs. 36%), GpG 3 (45% vs. 19%), and 
GpG 4 (30% vs. 15%) (Fig. 2). 

Increased VHL promoter methylation in malignant tumors 
and PGL 

The results from VHL methylation analysis were compared 
with tumor and patient characteristics (Table SI). This revealed 
higher Metl values in malignant as compared with benign 
tumors (mean 63% vs. 55%; P = 0.009). Furthermore, we found 
higher Metl in PGLs vs. PGGs (mean 66% vs. 55%; P = 0.001). 
No significant correlations were found between Met! values and 
gender (P = 0.5), age at surgery (P = 0.65), tumor weight or size 
(P= 0.4), epinephrine levels (P= 0.14), norepinephrine levels (P = 
0.53), dopamine levels (P = 0.7), constitutional SDHB mutations 
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Table 2. Results and correlations of methylation analyses of the VHL promotes 



Mean methylation 


Parameter 


CpG 1 


CpG 2 


CpG 3 


CpG 4 


Metl 


Tumors vs. adrenals 


Tumors 


79.9 


72.1 


44.2 


30.0 


56.6 


Adrenals 


37.3 


36.4 


18.9 


15.3 


27.0 


P value 


<0.0001* 


<0.0001* 


<0.0001* 


<0.0001* 


<0.0001* 


PGLs vs. PCCs 


PGLs 


88.7 


81.7 


48.9 


45.7 


66.3 


PCCs 


78.6 


70.6 


43.5 


27.5 


55.0 


P value 


0.002* 


0.003* 


0.293 


0.001* 


0.001* 


Malignant vs. benign 


Malignant 


85.6 


79.3 


46.7 


38.9 


62.7 


Benign 


78.6 


70.6 


43.7 


27.9 


55.2 


P value 


0.097 


0.011* 


0.482 


0.004* 


0.009* 


Methylation vs. gene expression 


l/H/. (assay 1)R = 


-0.35 


-0.31 


-0.24 


-0.26 


-0.36 


P value 


<0.0001* 


0.0003* 


0.0071* 


0.0027* 


<0.0001* 


Methylation vs. gene expression 


VHL (assay 2) R = 


-0.23 


-0.33 


-0.10 


-0.24 


-0.27 


P value 


0.009* 


0.0001* 


0.252 


0.006* 


0.002* 



*Statistically significant (P < 0.05); R, Pearson correlation coefficient. 



(P=0.l) or clinical MEN 2A syndrome (P= 0.31). CpG 4 (Chr3: 
10183220) showed the largest difference in methylation levels 
between PGLs and PCCs (mean methylation level 46% vs. 28%, 
P < 0.001) and between benign and malignant tumors (39% vs. 
28%,/'= 0.004) (Table 2). 

Suppression of VHL gene expression 

We studied the VHL gene expression in 74 tumors together 
with 13 normal adrenals using two individual VHL assays with 
qRT-PCR. To identify endogenous controls, we assessed 32 
candidate housekeeping genes. We found 5 genes to be stably 
expressed among PCC, PGL and normal adrenal references. 
These include: B2M, CASC3, MRPL19, PSMC4, and UBC 
(Fig. S2). Significant changes were found in tumors as compared 
with normal adrenals. In general, the VHL gene expression was 
lower in tumors than in normal adrenals (P < 0.0001 for assay 
1 and P = 0.002 for assay 2, Fig. 3). Furthermore, we found 
significantly lower VHL gene expression on PGLs as compared 
with PCCs for assay 2 (P = 0.047) but not for assay 1 (P = 0.72). 
Regarding tumor weight, we found significantly lower VHL 
expression in smaller tumors (Pearson correlation test R = 0.24, 
P = 0.007 with assay 2, R = 0.11, P = 0.2 with assay 1). Elevated 
norepinephrine levels were found significantly correlated with 
higher VHL gene expression (P = 0.02 for assay I, P = 0.15 for 
assay 2). No significant correlations were found between VHL 
gene expression and gender (P = 0.8 for assay 1, P = 0.5 for assay 
2), age at surgery (P = 0.58 for assay I, P = 0.73 for assay 2), epi- 
nephrine levels (P = 0.5 for both assay 1 and 2), dopamine levels 
(P = 0.9 for assay I, P = 0.3 for assay 2), constitutional SDHB 



mutations (P = 0.1), or clinical MEN 2A syndrome (P = 0.8 for 
assay I, P = 0.6 for assay 2). Furthermore, no significant correla- 
tion was found between VHL gene expression and malignancy 
(P = 0.97 for assay 1,P= 0.63 for assay 2). 

Inverse correlation between VHL methylation status and 
mRNA expression 

To further investigate the effect of VHL promoter methylation 
on its gene expression we compared Metis for the promoter with 
VHL mRNA expression levels as determined by each of the two 
qRT-PCR assays. Pearson correlation analysis showed significant 
correlation (R = 0.796, P < 0.0001) between the two mRNA 
expression assays supporting that no methodological errors are 
present since the same gene was analyzed with dissimilar assays. 
Furthermore, Metl values were inversely correlated with VHL 
expression (Pearson correlation test R = -0.359, P < 0.0001 for 
assay 1 and R = -0.232, P = 0.0077 for assay 2). With regard to 
the individual CpGs, all sites showed inverse correlation between 
methylation level and VHL expression for assay 1 (Table 2). 

Discussion 

Using a quantitative approach we assessed the methylation 
status of the promoters of 11 PCC/PGL susceptibility genes in a 
large cohort of PCCs, PGLs and adrenal reference samples. We 
found significant hypermethylation of the VHL gene promoter 
in the tumor samples compared with normal adrenal references. 
VHL gene expression correlated inversely with the epigenetic 
status of the tumors, showing lower expression in tumors and 
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Figure 1 . Comparison of VHL promoter methylation In normal adrenals and benign or malignant tumors. (A) Box plot showing Methylation index (Metl) 
for VHL in the different sample groups. (B) Column chart showing individual Metl for each sample. *WHO malignant cases. 



higher expression in adrenal references. Given that gene silenc- 
ing by promoter CpG methylation is an important epigenetic 
mechanism in tumorigenesis,'^ our results indicate that VHL 
promoter hypermethylation could play a role in the development 
of PCC and PGL. Welander et al. have shown lower expression 
of the VHL gene in tumors with loss of one gene copy." Other 
mechanisms for VHL inactivation could encompass involvement 
of microRNAs,"* mutations or, as our results indicate, promoter 
hypermethylation. Taken together our findings suggest that epi- 
genetic VHL silencing is a frequently occurring and potentially 
disease-related event in PCC and PGL development. 

In this study we observed very low methylation levels for all 
tested susceptibility genes except VHL. This is essentially in 
agreement with previous, non-quantitative publications discuss- 
ing methylation analysis of PCC/PGL susceptibility genes. "'" '^ 
However, our observations contrast previously published con- 
clusions concerning the VHL gene. This discrepancy may have 
several plausible explanations. Most important are the method- 
ological differences; the previous studies that report unmeth- 
ylated VHL gene promoters employed MSP"'"'"' and analyzed 
sequences located 127 nucleotides (6 CpG sites) downstream 
of the sequence analyzed in this study (Fig. 4). We found that 
the methylation density varies within the sequence of 4 CpG 



dinucleotides covered by our VHL assay (Table 2; Fig. SI). 
It is therefore possible that the previously assessed CpG sites, 
located downstream of our assay, are indeed unmethylated — as 
surmised by Welander et al. and Dammann et al. However, a 
methylated VHL promoter has previously been found by Moore 
et al. in renal cell carcinomas employing primers targeting a 
sequence 302 nucleotides (38 CpG sites) downstream of the 
sequence analyzed in our study (Fig. 4)." Previously, a genome- 
wide methylation analysis (Illumina Infinium HM27 DNA 
methylation assay) was undertaken for 145 PCCs and PGLs'^ 
which did not report hypermethylation in the VHL promoter. 
This discrepancy with our data most likely stems from the fact 
that the CpG sites covered by the Illumina 27K chip are dissimi- 
lar to those analyzed in this study. The nearest location for any 
individual CpG site covered by the Illumina assay is found -1 kb 
downstream of the sequence assessed herein. 

The increased VHL Metl levels in the tumors compared 
with normal adrenal references suggest that his is an important 
event in tumor development. A microarray study comparing rat 
PCCs with normal adrenal medullas revealed a number of dif- 
ferentially expressed genes, but did not identify VHL as one of 
them, indicating a difference between rat and human PCCs.^° 
We also observed significantly increased Metl levels for VHL 
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in malignant as compared with benign 
cases and the VHL gene expression was 
inversely correlated with Metl values. 
However, the VHL gene expression 
did not show any statistically signifi- 
cant difference between malignant 
and benign cases. Likely, the varia- 
tions in VHL gene expression may be 
present but can be too small to detect. 
Previously, gene expression profil- 
ing studies in human have not identi- 
fied VHL as a differentially expressed 
gene between malignant and benign 
PCCs,^''^'' which is in agreement with 
our observations. 

A CpG hypermethylator phenotype 
has been described for the PGL and 
5£)//5i2,i3 ^5 ^5 SDHx'* mutated 

entity. In this context it is of interest 
to note that the two cases with RET 
hypermethylation were malignant PGL 
harboring SDHB mutation (PGL-102 
and PGL-113). Neither VHL promoter 
methylation nor VHL gene expression 
correlated with SDHB mutations and 
accompanying CpG island methylator 
phenotype (CIMP).'' The possible 
association to clinical VHL syndrome 
could not be evaluated as this was only 
observed in one patient. 

In summary, our results show that 
the VHL gene promoter is hyper- 
methylated in pheochromocytomas 
and paragangliomas, accompanied by 
stunted gene expression, a finding that 
highlights the importance of VHL inac- 
tivation in the etiology of these tumors. 

Materials and Methods 



Tumor samples and controls 

Fresh frozen tissue material was obtained from the Biobank of 
the Karolinska University Hospital, Stockholm, Sweden. These 
samples had been collected with written or oral informed consent, 
as documented in patients' files, and sampling was approved by 
the local ethical committee. All tissues included were dissected 
and handled in accordance with established biobank protocols. 
The cohort included 96 tumors (83 PCCs, 13 PGLs; representa- 
tive for the ratio of PCCs vs. abdominal PGLs operated at the 
Karolinska University Hospital) (Table SI). Thirty samples of 
non-tumorous adrenal tissues were used as normal adrenal refer- 
ences (Table S2). These samples were obtained in conjunction 
with renal or adrenal surgery and histopathologically verified 
as non-tumorous by two expert pathologists independently. In 
addition, we included 4 samples of commercially available DNA 
from normal adrenal medullas (NAM-1— 4, Table S2) obtained 




CpG 2 CpG 3 CpG 4 Metl 
^ Adrenals # Benign ▲ Malignant 



Figure 2. Methylation at each of the 4 individual CpGs of the VHL promoter in normal adrenals and 
tumors. The scatterplot shows Pyrosequencing values for each CpG and case analyzed together with 
the Methylation index (Metl). Normal adrenals, benign and malignant tumors are indicated with col- 
ors according to legend. 



from autopsy of patients deceased from unrelated causes 
(Clinomics Biosciencies Inc.). 

Cases were characterized for gender, age of onset, hormonal 
status, existence of multifocal tumors, occurrence of multiple 
or bilateral tumors, tumor weight and tumor size, recurrence, 
duration of follow-up, patient outcome and time for metasta- 
sis. In addition, occurrence of constitutional SDHB mutations 
and clinical VHL, MEN 2A, and NFl syndrome is given in 
Table SL In the current paper malignant cases were identified 
based on metastasis and/or extensive local invasion according to 
the Armed Forces Institute of Pathology (AFIP) guidelines,^'' due 
to higher sensitivity in detecting potentially deleterious tumors 
(herein referred to as malignant). Additionally, we have given 
the classification according to the World Health Organization 
(WHO) guidelines' (referred to as WHO malignant), since these 
are also widely accepted for the classification of these tumors. The 
WHO criteria base malignancy solely on the presence of distant 
metastases.' A small representative piece was cut from 115/130 
(88%) of all analyzed samples. The piece was thawed and fixed 
in formaldehyde solution, embedded in paraffin, sectioned and 
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stained with Haematoxylin and 
Eosin. Tumor or normal cell content 
was subsequently determined by two 
expert pathologists, independently 
revealing representative tumor or nor- 
mal cell content of more than 70% in 
92% of the included cases (data not 
shown). Cases that revealed tumor 
or normal cell content <50% were 
excluded from the study. 

Preparation of DNA and RNA 

Extraction of genomic DNA was 
performed using DNeasy Blood and 
Tissue DNA isolation kit (QIAGEN 
AB). Quantification and initial 
quality check of the extracted DNA 
was performed using a NanoDrop 
AlOO (Thermo Scientific). EpiTect 
Bisulphite conversion kit was used 
for bisulphite conversion of 500 ng of 
genomic DNA from all tumors and 
normal samples (QIAGEN AB). The 
method is described in detail in a pre- 
ceding publication.^"^ 

Extraction of total RNA was 
performed using the wzVVana 
miRNA Isolation Kit (Ambion, Life 
Technologies, Applied Biosystems) 
with no changes to the original proto- 
col. RNA quality was measured using 
an Agilent 2100 Bioanalyzer (Agilent) 
with the RNA Nano 6000 chips and 
reagents (Agilent) (Tables SI and S2). 

Pyrosequencing and data analysis 

Pyrosequencing was performed to 
assess the methylation level of PCC 
and PGL susceptibility genes in 96 
tumors (83 PCCs, 13 PGLs) and 
34 adrenal /adrenal medullary refer- 
ence samples. The assays for EGLNl, 
MAX, MENl, NFl, RET, SDHAF2, 
TMEM127, and VHL are commer- 
cially available (PyroMark Assay 
Database, Geneglobe, QIAGEN). 
The SDHB, SDHC and SDHD 
assays were ad hoc-designed using 
PyroMark Assay Design software 
version 2.0 (QIAGEN). Both newly 
designed and commercially available 
assays were extensively tested for assay 
bias as previously described'^ (results 

not shown). Assay specifications and sequences are summarized 
in Table S3. We also made several attempts at designing primers 
for SDHA, and for the VHL region corresponding to the target 
sites in an article by Dammann et al.,"" however these assays did 
not produce reliable results and were therefore excluded. 
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Figure 3. Comparison of promoter methylation density and gene expression for VHL. (A) Boxplots show- 
ing relative VHL mRNA expression (y-axis) In the sample groups (x-axIs) determined by qRT-PCR using 
two different assays. (B) Scatterplot Including trend line showing Metl values (x-axIs) and relative mRNA 
expression levels determined by assay 1 (y-axIs). 



Pyrosequencing has been described in detail in a previous 
publication.^^ In short, pre-made master mix — including DNA 
polymerase from QIAGEN — was used to amplify 1 ^jlI of bisul- 
phite—treated DNA in a PGR reaction, using the following con- 
ditions: 95 °G for 15 min, followed by 45 cycles of: (95 °C 30 s, 
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3p25.3 
VHL 

ATG 

Exon 1^ ^ 

Nucleotide number 1 841 1180 1181 1303 1304 1484 

I I 
I I 

" « 

^^"^ Bisulphite treated, methylated sequence """"""--.^^ 

TGAGTGTTTATGTTTGTAGTTTTAATTGTTTTGAAGTAAATATTTGAI I 1 1 1 I AATTTTTATTAGAGTG 

TTTTGTATATAGTAGGTTTAATTAI I I 1 1 I 1 1 I 1 1 1 lATTAATTATTTATGTTTTGTAA GAATTTAGTTA 

GTTGAI I I 1 1 IGTATmATAAGCSTGATGATTGSGTSTTTTCGTGTGASATGCGTTATTTTCGAATTT 

TGTTACGACGTCGGTATATTGCGCGTTTGATATGAAGAAAAAAAAAATTTAGTTAGTTTATTAGGTA 

TAGTGGTTAAGGmGTAATTTTTGTATTTTGAGAGGTTAAGGTAGGAGGATTATTTGAATTTAGGA 

GTTCGAGATTAGTTTAGGTAATATAGCGAGATTTCGTTTTAAATAATAAATAAAAATAATTAGTCGG 

GTATGGTGGTGCGCGTTTATAGTATTAATTATTCGGGASSTTSAGGCGAGACGATCGTTTGAGTTA 

GGGAGGTTAAGGTTGTAGTGAGTTAAGT ATAGAGTGAGA 

TmGTTTTAAAAAAAAAAAAAATATTAAATTTTAGAGGGGCGAAAAAAAATTTTATAGTGGAAATA 

TAGTAACGAGTTGGTTTAGTTTCGTTTTCGTTATAACGGTTTACGGTGT 

STATAGTTTTCgGTCSGTTATTTTCGCGAGCGCGTTTTAI I 1 1 I lATCGAGCSCSCGCGAAGATT 
ACGGAGGTCGATTCGGGAGCGCGTACGTAGTTTCGTTTCGCGTTCGATTCGCGGAT 

GGGTGGTTTGGATCGLL .GTTTCGGAGGGCGGAGAATTGGGACGAGGTCGA 
GGTAGGCGCGGAGGAGGTAGGCGTCGAAGAGTA CGGTTTTGAAGAAGACGGCGGGGAGGAGTCG 
GGCGTCGAGGAGTTCGGTTCGGAAGAGTTCGSTTCGGA66AATT666CGTC6AGGA66AGATGG 
AGGTCGGGCGGTCGCGGTTCGTGTTGCGTTCGGTGAATTCGCeCGAGI 1 1 1 1 1 lAGGTTAI I 1 1 1 IG 



Figure 4. Schematic overview of the VHL promoter region. Underlined 
letters in black represent the annotated promoter sequence (GenBank 
accession No AF010238), letters in green denote the sequence ana- 
lyzed in this study, letters in yellow denote the primer sites analyzed by 
Dammann et al.,"" letters in blue represent the sequence analyzed by 
Moore et al.," bold letters represent CpG islands and letters in red mark 
the translation start site. 



56 °C [all assays except SDHB, SDHC, and SDHD; 53 °C] 30 s, 
72 °C 30 s), and 72 °C for 10 min. The PGR product was then 
run in a PyroMark Q96 ID pyrosequencer, employing unique 
run files for each separate assay (available for download through 
QIAGEN). Quantification of methylation density was done 
using the PyroMark Q96 software. Negative controls (water sam- 
ple, no primer sample, no DNA sample and 0% methylated com- 
mercially available DNA [QIAGEN]) and positive controls (50% 
and 100% methylated commercially available DNA [QIAGEN]) 
were included in all experiments. 

For each assay and sample, a methylation index (Metl) was 
calculated as the mean methylation of all GpG dinucleotides 
covered in the assay. Subsets of cases were re-analyzed with 
repeated PGRs and subsequent Pyrosequencing to validate the 
initial results. These tests had similar end results (±2%) for all 
GpG dinucleotides in all cases (data not shown). Subsets of cases 
(including both PGGs and PGLs) were bisulphite treated a second 
time and subsequently analyzed using PGR and Pyrosequencing, 
with analogous end results (data not shown). DNA from mul- 
tiple cases was also re-extracted for the same experiments and 
Pyrosequencing showed the same end results (±2%) for all GpGs 
in all cases tested (data not shown) , indicating the robustness and 
reproducibility of the methods employed. 

Quantitative reverse transcription real-time PGR 
(qRT-PCR) 

qRT-PGR analysis was performed for 74 of the tumors, 
together with 13 of the non-tumorous adrenal reference samples 



where RNA was available. The 4 normal medulla samples were 
not included. cDNA was synthesized from 100 ng total RNA 
using High Gapacity cDNA Reverse Transcription Kit, fol- 
lowing manufacturer's protocol (Applied Biosystems). For the 
qRT-PGR, 1 |jl1 of cDNA template was used together with PGR 
reaction mix and the particular TaqMan gene expression assay: 
VHL mRNA assay 1 (Hs00184451_ml, Applied Biosystems) 
and VHL mRNA assay 2 (Hs01650959_ml) targeting total 
VHL transcripts. 

As recommended by the Minimum Information for 
Publication of Quantitative Real-Time PGR Experiments 
(MIQE) guidelines,^** we tested 32 different endogenous con- 
trol genes in triplicates for 6 different PGG, PGL, and normal 
adrenal tissue samples. These genes are included in the TaqMan 
Array Human Endogenous control plate (Applied Biosystems). 
Results from the testing of the endogenous control run are sum- 
marized in Fig. 82, and their specifications are summarized 
in Table S3. The 5 most stable endogenous controls chosen 
for further experiments were selected from the ones displaying 
the smallest variation among the samples i.e., B2M, CASC3, 
MRPL19, PSMC4, and UBC. The geometrical mean of the 
expression of these genes were used to normalize all expression 
data. Triplicates of PGG, PGL, and normal adrenal samples 
were run together with endogenous controls and non-template 
controls in 384-well plates using standard quantitative method 
of amplification in a 7900HT Real-Time PGR System (Applied 
Biosystems). Relative expression levels for the VHL gene were 
reported as 2"^'-'^. 

Statistics 

Pearson correlation test was used to compare methylation 
density and gene expression for methylated genes and to compare 
age at surgery with methylation density and gene expression. 
Unpaired student f-test was performed to analyze differences in 
the methylation and gene expression levels between the tumor 
groups and normal adrenals, respectively. Unpaired student 
^-test was also used to compare gender with methylation status 
and gene expression. Gomparisons between tumor subgroups 
were always done employing the AFIP criteria for malignancy. 
When relevant, WHO malignant cases are separately specified. 

A P value of < 0.05 was regarded as statistically significant. 
The statistical analyses were performed using STATISTIGA vlO 
(STATSOFT) and Figures 2 and 3B were generated with Think 
Gell (think-cell Software GmBH) 
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